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Adiabatic quantum computation is based on the adiabatic evolution of quantum systems. We 
analyse a particular class of qauntum adiabatic evolutions where either the initial or final Hamilto- 
nian is a one-dimensional projector Hamiltonian on the corresponding ground state. The minimum 
energy gap which governs the time required for a successful evolution is shown to be proportional to 
the overlap of the ground states of the initial and final Hamiltonians. We show that such evolutions 
exhibit a rapid crossover as the ground state changes abruptly near the transition point where the 
energy gap is minimum. Furthermore, a faster evolution can be obtained by performing a partial 
adiabatic evolution within a narrow interval around the transition point. These results generalize 
and quantify earlier works. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

Quantum adiabatic evolution starts with the ground 
state \s) of the initial Hamiltonian Hj in an A^- 
dimensional Hilbert space, and evolves it slowly enough 
to the ground state \t) of the final Hamiltonian Hf. The 
evolution uses the time-dependent Hamiltonian 

= (l-/i)H, +AiHt , fie [0,1]. (1) 

The parameter /i is a function of the time r. The quan- 
tum adiabatic theorem p| bounds the total evolution 
time r required for a successful evolution. Let the eigen- 
spectrum and the excitation gap of be 

Eo,fi < Ei,fj^ < • ■ ■ < EN~i,fj. , 

9p. = Ei^f^ - i?o,Ai • (2) 

The adiabatic theorem states that one can reach \t) with 
probability close to 1, when 

r > e {g^^J\ - Htll) , 5min = min^ g^, . (3) 

Conventionally, the Hamiltonians are normalized such 
that II — Htll — 9(1), and F is bounded from below 
essentially by g^f^- Thus the knowledge of the minimum 
energy gap gmXn is essential to determine the minimum 
time for successful evolution. In general, estimating gmin 
is not an easy task but it can be estimated for some spe- 
cial cases BMlSSB- 

In this paper, we analyze a special case when Ht = 
— |t)(t| is a one-dimensional projector Hamiltonian on 
its ground state \t). Such kind of projector Hamiltoni- 
ans naturally appears in solutions to decision problems. 
In Section II, we analyse the eigenspectrum of H^ with 
Ht = — |t)(i|. We show, under certain assumptions re- 
garding the eigenspectrum of H^, that gmin scales as the 



overlap a — \ {s\t) \ of the ground states of Hs and Ht. We 
also derive the expression for /Zmin where = gmin- For 
5min = 0(a), ([3]) implies that F = 17(a~^). In Section HI, 
we show that the ground state of H^ evolves significantly 
only within a narrow interval [/i^, /i^] around /Xmin- Elx- 
ploiting this property, we present a partial adiabatic evo- 
lution algorithm with the time complexity F' = Q,{a~^), 
which is faster than the standard adiabatic evolution. In 
Section IV, we conclude by discussing the relation of our 
work to earlier works on this subject. 

II. MINIMUM ENERGY GAP 

To calculate the minimum energy gap, we first find the 
eigenspectrum of H^. With Hf — — |t)(t| in ([T]), we have 

H^ = {l-^J)Hs-^l\t){t\. (4) 

The eigenspectrum of above Hamiltonian can be analyzed 
in a similar way as the eigenspectrum of corresponding 
unitary operator was analyzed in We work in the 
eigenbasis of Hs, chosen such that 

(^|Hs|£)=e^ , Q = io<ii<---<S,N-i ■ (5) 

For simplicity, we consider \s) = \£ — 0) to be the non- 
degenerate ground state of Hs. We make the following 
assumptions regarding the eigenspectrmn of Hsi 

\{.3\t)\ = a « a , it 1 , IIHsll = ^N-i > 1. (6) 

The first one can always be satisfied by appropriately 
scaling Hg (and hence a). The time needed to distinguish 
the ground state of Hg from the excited states, 17(1/^), 
is then much smaller than the time scale of the algorithm, 
F = Q{l/a). The second one constrains ||Hs|| relative to 
the initial excitation gap ^ , and the third one constrains 
IIH.II- 

Let |i?fe,;i) be the normalized eigenvectors of H^ with 
eigenvalues E^^^. We have 
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H^l^fc,^) = Ek,t,\Ek,^,) - [(1 - /i)Hs - M|i)(i|]|Sfe,^). 
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Left multiplication by (£| and (£|Hs = gives 

= ^i){^\Ek^^) - ^Ji{^\t){t\Ek,^:). 



Thus 



k,fi) 



m{t\Ek., 



(7) 



It gives 



= ^^{t\Ek,^)J2■ 



\m\ 



6(1 - ^J.) ~ Ek^i^ 
and we find the secular equation for to be 



E 



6(l-M)-i?fc„ 



1 



(8) 



(9) 



Since 6 ^ 0, the L.H.S. of above equation decreases 
monotonically as Ek,fj, decreases from to — oo. On the 
other hand, the R.H.S. is fixed, so the equation can have 
at most one negative solution for Ek^fj,. We will see that 
above equation has a unique negative solution, which is 
obviously the ground state energy Eq^^ of H^. 

We assume that the two lowest solutions of © obey 



|i?fc, J « (1 - /i)6 



(10) 



To find them, we Taylor expand the i ^ contribution 
in ^ and ignore 0{E'^ ^) terms. That results in the 
quadratic equation, 



2 77-1 



A,. 



BlEu., 







(11) 



yielding two solutions consistent with l-E/t.^l ^ (1 — /i)6 



The coefficients i?^ are 



A„ 



where 



Ti 



E 

^#0 



- - , B„, = 



/To 



\mi 



pe{i,2}. 



(12) 



(13) 



We note the bounds Cw^i 1^ ^ arising from 

Y.iW\t)? = 1- Also, putting xi = \{l\t)\ and ye = 
1(^1^)1/6 the Cauchy-Schwartz inequahty i^xiyiY < 
i2xeT>yh we get 



T?<E 



< T2 



(14) 



The two solutions i?±,p of (fTT|) have the product 



= -a'/Bl. Hence 



E±,^ = ± — (tan??^) 



±1 



±^^(tan,,)±^ 



(15) 



The sum of the two roots determines the angle 1]. We 
have E+^^ + E_^^ = -A^/Bl = -(2a/B,J cot 27?^. Thus 

A,, 



cot 27]^ = 



1 



2aB^ 2aVT2 \ M 



Ti- 



1-M 



(16) 



with 77^ e [Oif]- As 77^ is positive, i?-,^ is indeed the 
unique negative solution of ([9]) and hence the ground 
state energy Eq^^ of H^, while E+^^ is the first excited 
state energy of H^. 

With Ti > 0, let us define the crossover point fi* , and 
deviation from it e as 



1-M* 



= Ti 



1 



M = 



£ = 1 - 



Z^" 1 + Ti fi 

By definition, A^* — and r/^* = j. We also have 

1 / 1 - /i* 1 - £ 



(17) 



cot 2?/^ = 



2aVT^ V 
(1 + Ti) ^^ 
2aVTl ^' 



1 



(18) 



The bound T2 < ^ and the assumption a ^ 6 give 
aVT2 < a/Ci *C 1- Then | cot 277^ | is large for e not close 
to 0. On the other hand, for ^ close to /i*. 



|e| <C 1 : cot 277^ = Lu{^ - fi*) 



(1 + Tr)! ^ a ^ 



From p5|) . we obtain the excitation gap as 
a{l - fj.) 



(19) 



5p 



(tan 7/^ + cot?7p) 



2a(l - n) 



CSC 277^ . 



(20) 



Since esc 27;^ > |cot27;^| ^ 1 for /i not close to /i*, 
is close to its minimum only when /i is sufficiently close 
to /i*. The size of this region is characterized by the 
parameter lo. Explicitly, using (jl9p for 11 ~ 11* <^ 1, we 
get 

2a(l - - (/^ - /^*)) 77- — 27 



2a(l-/i*) 



1 - /i* 2 ^ V 



At its minimum. 



2a(l -^*) 



1 - 



1 



2Lu^{i~^i*y 



(21) 
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With u! ^ 1, the deviations of these values from their 
values at the crossover point are tiny. The assumption 
^ 1 gives \/T^ > ^^Li = ^^(1), and hence 
= 0{a). Also, for ^ close to fi*, and (^1]) 

can be combined as — gmin esc 277^ . 

As ^niin — 0(a), ([3]) implies that the time required 
for the standard adiabatic evolution to be successful is 
r > 0(a~^). We observe that the state \t) can also be 
obtained by a simple scheme of 0(a~^) times prepara- 
tion and subsequent measurements of the state |s) in a 
suitable basis. Hence, the standard adiabatic evolution 
does not give any speedup over the simple scheme. In 
the next section, we show that if we know the crossover 
point fi* then we can achieve a faster algorithm with the 
time complexity 0{a~^). 



[/i",/!"*"]. We will also find below that the ground state 
of changes substantially only within the inter- 
val This property can be used to construct 
a faster adiabatic algorithm which performs a partial 
adibatic evolution only within the interval and 
safely skips the evolution outside this interval. 

The validity of our analysis relies on PH)) . With = 
li?-!-.,,,! -I- \E-_„] and ((20)) . this vahdity condition holds 
< Ci\/TI/2. Now for n e 

and > ii/2^N-i ^ 1 



provided acsc 277^ 
a CSC 277^ < aVT+c? 



due to the assumption ([6]). Since a <C ^i, the validity 
condition can be satisfied by keeping c small compared 
to £,i/a. That also keeps the interval narrow, 
with c <C w as per (jl9p . Thus our analysis is valid within 
the interval 
Evolution of 



is obtained by ((23 and (gS]), 



I.e. 



III. PARTIAL ADIABATIC EVOLUTION 

Before presenting the faster algorithm, we first com- 
pute the overlap of the ground state E-^^ of with 
the initial and final ground states, \s) and \t). With the 
normalization condition \ {£\E^^^)\'^ = 1,0 gives 

12 1 



E 



(22) 



For E^^^ < - m), 0(£;-,p/CHl - ^^)) terms can be 
ignored in above equation to get 



\{s\E-,i,) \ = sin77p , 
|i?-,M>l = [(i-l)/x/T;]cosr;^ 



(27) 



We have 77^ w ^ for < fi^ , and ry^ w for /i > So 
in passing through the interval [fi" , fi'^] , the ground state 
transforms from being very close to the initial 
state \s) to being almost orthogonal to |s). (Note that 
when rji^ is a smooth function, \E-^^) is close to \s) for 
all < /J,~ , even though our analysis does not hold for 
all fi < fi^ ■) Simultaneously, the overlap of with 
the target state \t) increases from zero to 



Ei 



E 



l(^|i>l' 



1 



(23) 



{t\E_ 




(28) 



or 



E' 



\E-,,W 



(24) 



where we have used the definition of . Using ([T 

for E^^fj^, above equation gives 



COS^yy^ /i2|(t|i?_,^)|2 



mE-.,)\ = 



1 COS 7?^ 



(25) 



\{t\E_.^)\ = U-\)/^2]coBi^^ • 



To compute |(s|i?_,p)|, we put ^ = = s, ^0 = 0, ([T5 
and in © to get 



l(s|i?~,M>l=/^- 



sm 



(26) 



Now, consider the narrow interval where 
/i^ = ^* ± cw"^ with 1 ^ c ^ a;. Outside this in- 
terval, we have |^ — > ccj^^, and | cot2?7p| > c from 
(fTO)) . Therefore, r]^>^+ is close to zero, 7y^<^- is close to 
■| , and (7^ = (7inin CSC 2/]^ is much larger than gmin outside 



We have 1 > Ti/a/T2 > £,i/^n-i due to the bounds 
mentioned after (fT3|) . The assumption ^i/^jv-i ^ 1 then 
implies that \t) has a significant overlap with \E_^^+). 
Hence \t) can be obtained by few preparations and sub- 
sequent measurements of \E_^^+). 

We now define the partial adiabatic evolution as evo- 
lution from H^- to H^+, as opposed to the complete adi- 
abatic evolution from Hg to Ht. The resultant algorithm 
executes the three steps below: 

(1) The initial state \s) is prepared as the stable ground 
state of the Hamiltonian H^. At r = 0, the Hamiltonian 
is suddenly changed to H^- , without disturbing the state 
Is) [J. The system is then in the ground state 



with probability sin^ ' 



(2) The Hamiltonian evolves from H^- to H^+, linearly 
in time over duration T. The system encounters the min- 
imum excitation gap gmin during this evolution, and the 

with probability 



reaches the state 



E_ 



ymin 



state 

close to 1 for T > 2cuj~ 
(3) The state of the system is measured. The state 
\E_^^+) yields the target state \t) with probability 
TI/T2- These three steps are repeated until we find \t). 

The combined success probability of the three steps 
is Pad = sin^ 77„- T2/T2. The overall complexity of the 
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algorithm is, to leading order, 

(29) 

making use of and (PT|) . Since Tp e [S.N-n^i ^^'i 
we have assumed both and £,n-i to be 6(1), the fac- 
tor tI''^/T\ is also 9(1). That makes V > &{a~^), and 
hence the partial adiabatic evolution is indeed quadrati- 
cally faster than the 0(a^^) classical search algorithms. 

IV. DISCUSSION 

We can obtain Roland and Cerf 's results [7] as a special 
case of our partial adiabatic algorithm. There Hj — In — 
\u){u\, \u) = J2j \j)/VN and a = 1/VN. We then have 
e,^, = 1 and Ti - T2 = Y^i^^ = I - c? . It 

follows that the crossover point is fi* ~ 1/2, with w k, 
2/ a and g-min ~ ol. The width of the narrow interval 
is 2ca;^^ = ca ^ 1 as desired. Using the partial 
adiabatic evolution algorithm, we obtain the target state 



in time V > cTj^^/aTf « c\/N, which is optimal up to 
a constant factor. Roland and Cerf obtained the optimal 
algorithm by performing local adiabatic evolution which 
performs the evolution slowly around the crossover point 
where the energy gap is small. Note that the knowledge 
of the crossover point /^* is essential to get the optimal 
algorithm, either by partial evolution or local evolution. 
Another special case of our analysis is due to Farhi et 
al. [2], where Hs is a sum of single qubit Hamiltonians. 

By time reversal symmetry, our analysis can be ex- 
tended to the problem where Hg = — |s)(s| and Hj is 
a general Hamiltonian. The required interchanges are 
t ^ s and ^ 1 - M. Also, Tp - IO»lV?^ 
where represent the eigenspectrum of Hj with 

^j^t = 0. The crossover point becomes fj,* = Ti/(1-|-Ti), 
and {g^, uj, T} can be calculated. Znidaric and Horvat Q 
have studied a particular case of this type, with Hj rep- 
resenting instances of an NP-complete problem. 

Acknovifledgments: I thank Prof. Apoorva Patel for 
going through the manuscript and for useful comments 
and discussions. 



[1] A. Messiah, Quantum Mechanics, Chapter XVII (Dover, 

New York, 1999). 
[2] E. Farhi, J. Goldstone, S. Gutman, and M. Sipser, e-print 

'quant-ph/0001106 
[3] E. Farhi, ,]. Goldstone, S. Gutmann, and D. Nagaj, e-print 

quant-ph/0512159 
[4] L. M. loannou and M. Mosca, International Journal of 

Quantum Information, Vol. 6, No. 3 (June 2008). 



[5] Z. Wei and M. Ying, e-print [quant- ph/06 040'f7l 
[6] M. Znidaric and M. Horvat, Phys. Rev. A 73, 022329 
(2006). 

[7] J. Roland and N. J. Cerf, Phys. Rev. A 65, 042308 (2002). 
[8] A. Tulsi, General framework for quantum search algo- 
rithms, e-print quant-ph/0806.1257. 



